Neutron diffraction study of TbMnOs: Magnetic structure revisited 
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Magnetic ordering in TbMnOs was revisited by a neutron diffraction study. In addition to the pre- 
viously reported j4-type and G-type modulated structures with a propagation vector (0, gMn,0)r'" 
we found C-type and _F-type orderings of the Mn moments. All the components appear below 
T^" = 46 K, and the G, G, and F-type components are enhanced below Tiock = 28 K, where 
^Mn ordering locks at its low temperature value. The locking of the propagation vector also yields 
the squaring up of the spin arrangement. The magnetic moments of Tb ions show a quasi long- 
range ordering below T^*^ = 7 K. It drastically promotes the development of the G, G, and _F-type 
components while suppressing the ^-type components. 



I PACS numbers: 75.47.Lx, 75.25.+Z 

Pcrovskitc manganese oxides iiMnOa {R = 
trivalent rare earth ion) are famous as the parent 
materials of the colossal magnetoresistive manganites. 

\ The Mn ions in i?Mn03 become trivalent having 
the ifgSg configuration. This electronic configuration 
induces a degree of freedom of the Cg orbital, which 
produces many novel properties. In the most well stud- 
ied material, LaMnOs, the Cg orbitals form staggered 
ordering of d^^^-r^ and d^y2_^-i orbitals, which induces 
the layered-type (so-called ^-type) antiferromagnetic 
ordering of Mn spins. In contrast, when the ionic radius 
of R is substantially small, /JMnOa shows a modulated 

' spin ordering with a propagation vector {Q,qMm^) (hi 
the orthorhombic Pbnm cell)^^*^*^ The emergence of 
the incommensurate spin structure is explained as a 
competition between the nearest-neighbor (NN) spin 
interactions and the ncxt-nearcst-neighbor (NNN) intcr- 

I actions caused by the combination of the GdFeOs-type 
distortion and the d^^'^ / d^y2 orbital ordering (We 

; will refer to this model as the NN-NNN model. )i^ 

I TbMnOa is one of the latter series of i?Mn03. It 
shows an incommensurate sinusoidal spin ordering with 
<?Mn ~ 0.295 below Tjv ~ 41 K with spins oriented along 

i the [010] direction, qun decreases as temperature is low- 
ered, then locked at qun = 0.28 below Tiock 30 Ki 
Very recently, it was found that the spontaneous elec- 
tric polarization P parallel to the c axis appears below 
Tiock Moreover, the magnitude or the direction of P 
can be drastically changed by applying a magnetic fieldi 
The magnetic control of the ferroelectric polarization will 
provide an attractive possibility of new magnetoelectric 
devices. 

These recent discoveries in the i?Mn03 series moti- 
vated us to reinvestigate the magnetic properties of these 



materials in detail. Here, we report the result of a neu- 
tron diffraction study on the elastic properties of the 
magnetic ordering in TbMnOa. The study of the inelas- 
tic properties by neutron scattering is underway, and will 
be pubhshed elsewherei^ 

A single crystal of TbMnOs was grown by the floating 
zone method. A detailed procedure of the sample prepa- 
ration was described elsewhereji The neutron diffraction 
experiments were performed using the triple axis spec- 
trometer GPTAS installed at the JRR-3M research re- 
actor in JAERI, Tokai, Japan. A neutron wave length 
of h = 2.57 A-i was selected by the 002 reflection of a 
pyrolytic graphite (PG) monochromator. The horizon- 
tal collimation of 40'-40'-40'-blank (from monochroma- 
tor to sample) was used for most of the measurements, 
while the 20'-20'-20'-blank collimation was utilized for 
the investigation of temperature dependence of the wave 
vector of the magnetic ordering. Two PG filters were 
placed before the monochromator and after the sample to 
suppress contaminations of higher-order harmonics. The 
crystal structure has the Pbnm symmetry, and the mea- 
surements were performed in the (Okl) scattering plane 
(a* = 1.08 and c* = 0.853 A'^). The sample was 
mounted in an aluminum can filled with helium gas, and 
was attached to the closed-cycle helium gas refrigerator. 
All the temperature dependence data were collected upon 
heating to avoid uncertainty due to the hysteresis. 

First, we surveyed magnetic reflections in the (Okl) 
scattering plane with the relaxed collimation (40'-40'-40'- 
blank). Figure ^ shows neutron scattering profiles along 
the (0, fc, 0) and (0, k, 1) lines at 4 K, 10 K, and 55 K. At 
4 K, surprisingly many superlattice reflections due to the 
magnetic moments of both Mn and Tb ions are observed 
at every position with k k, n/7 [n = integers). At 10 K, 
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FIG. 1: (Color online) Neutron difltraction profiles along (a) 
the (0, fc,0) line and (b) the (0, fc, 1) line at 4, 10, and 55 K. 
Vertical lines indicate the positions of the Mn (fc ~ integer ± 
2/7 and 6/7) and Tb (fc ~ integer±3/7 and 6/7) superlattice 
peaks. "Al" denotes diffractions from the aluminum sample 
cell. 



on the other hand, only the Mn peaks are observed with 
broad peaks coming from the short-range order of the Tb 
moments. 

The magnetic reflections due to the Mn ordering were 
observed at {h, k ± ^Mn, with qun = 0.28, where solid 
vertical Unes are drawn in Fig.Q] The observed Mn peaks 
can be classified into five groups depending on the values 
of h, fc, and I as ^-type {h + k = even and I = odd), 
G-type {h + k = odd and I = odd), C-type {h + k = odd 
and I = even), and f -type (h + k = even and / = even). 
Quezel et al. observed strong yl-typc reflections and 
weak G-type peaks, but they did not observe any C- 
typc and F-typc rcflcctionsii They concluded that the 
^-type peaks come from a sinusoidally modulated spin 
ordering with spins parallel to the b axis. In contrast 
to Quezel's result, we observed clear peaks of both G- 
type and F-type in addition to the ^-type and G-type 
peaks. Since the small number of the observed magnetic 
peaks prevents us to unambiguously determine the real 
space spin configuration, we just performed a prelimi- 
nary analysis of the magnetic structure: We calculated 
the intensities for the sinusoidal and the helical spin con- 
figurations assuming the spins lie on the ab, be, or ca 
planes, and compared them with the observed intensities 
at 10 K. The intensities of the ^-type peaks are best de- 
scribed by the sinusoidal model with spins parallel to the 
b axis as reported by Quezel et al. Analysis of the G-type 
peaks gave the similar result. For the G-type and i^-type 
peaks, however, we failed to satisfactorily reproduce the 
observed intensities, though the existence of the peaks 
on the (0, k, 0) line indicates the spin arrangements have 
substantial x or z components. 

Figures |2Ia)-(d) show temperature dependences of the 
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FIG. 2: (a)-(d) Temperature dependences of the intensities 
of the Mn superlattice peaks of (a) ^-type, (b) G-type, (c) 
-F-type, and (d) G-type. (e) Temperature dependence of 
the intensity of the Mn third harmonics at (0, 0.835, 1) = 

(0, 3(JMn, 1). 



scattering intensities of the Mn magnetic peaks. Each 
panel corresponds to the superlattice peak of (a) A-type: 
(0,0.278,1) = (0,gMn,l), (b) G-type: (0,0.721,1) = 
(0, 1 - QMn, 1), (c) f-type: (0, 0.280, 2) = (0, qun, 2), and 
(d) G-type: (0, 1.280, 2) = (0, 1 + qun, 2). All the inten- 
sities start to increase at T^" = 46 K, which is slightly 
higher than the reported transition temperatures 
There are two distinct anomalies in the temperature pro- 



files below T^": with decreasing temperature, the slopes 
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show upturns at Tiock = 28 K, and then a drastic drop of 
the ^-type peak (a) or increases of the others (b)-(d) be- 
low T'^^ ~ 7 K. The former temperature coincides with 
the temperature where the wave vector of the Mn order- 
ing locks at its low temperature value as described later. 
The latter anomaly is concomitant with the ordering of 
the Tb moments (sec below). All the weak (G, C, and F) 
components show a similar temperature dependence, and 
they gradually develop below Tj^". Their temperature 
dependence is quite different from the sharp increase of 
the ^-type component, suggesting that the weak compo- 
nents arise from some secondary effect. One of probable 
causes is the development of the ordering of the Tb mo- 
ments, which is manifested by the steep increase of these 
components in compensation for the ^-type component 
on the (quasi) long-range ordering of the Tb moments 
(Fig.©. 

In addition to the fundamental (Iqun) peaks described 
above, we also observed the third harmonics (SgMn) of 
each type of the Mn ordering at {h, k ± SgMn, (dashed 
lines in Fig. ^ . Though we could not detect the third 
harmonics for the F-type ordering, it may be because of 
the weakness of the fundamental (IgMn) F-type peaks. 
The existence of the higher harmonics of the Mn order- 
ing means that the spin ordering is not an ideal sinu- 
soidal wave. The ratio of the intensity of the 3gMn peak 
to that of the Iqym peak is ~10~^. This value is much 
smaller than that for an ideal square wave (1/9), indi- 
cating the deviation of the spin arrangement from the si- 
nusoidal wave is quite small. Nevertheless, the existence 
of the higher harmonics is interesting in that the NN- 
NNN model predicts the square ordering of the Ising-like 
spinsi^ In Fig.[2Ie) is shown the temperature dependence 
of the intensity of the SgMn peak of the ^-type ordering. 
It gradually develops below Tiock, indicating that the si- 
nusoidal ordering pattern of the Mn moments becomes 
distorted at T < Tiock- 

The superlattice peaks due to the ordering of the Tb 
moments were observed at {h,k ± gxbiO with h,k,l = 
integer and gxb = 0.42 (dotted-and-dashed lines in 
Fig. nj. The Tb peaks are observed in all zones simi- 
lar to the Mn ordering. Though the previous study on a 
single crystal reported the Tb peaks are broad and less 
intense compared to the Mn peaksji we observed quite 
sharp and much intense peaks at the lowest tempera- 
ture. Figure |2Ia) shows the temperature dependence of 
the Tb peak intensity at (0,0.424,1) = (0,gTb,l)- It 
shows a steep increase below T'P^ = 7 K. In addition, 
superlattice peaks arc observed at k ± 0.15 in all zones 
(dashed lines in Fig. ©. The intensity of these peaks 
shows the same temperature dependence as the Tb peaks 
[Fig-EJb)], indicating that these peaks are related to the 
Tb ordering. Since the positions of the superlattices cor- 
respond to 2gTb, these peaks may be attributed to the 
second harmonics of the Tb ordering. The existence of 
the even harmonics means a space symmetry breaking^ 
On the other hand, a magnetic ordering can induce the 
second harmonics of lattice origin^ The exact origin of 



o 






CD 




(a) : 


o 


4 


* • m n AOA 1 \ 


c 




\ =(0, qjb. 1) 


13 

o 


o 
O 


■ \ IT^" - 


o 




\ 


c. 


2 


It 


>. 






w 


1 


: \ ~ 


o 






— 







o 

CD 




1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

(b) ; 


o 


8 


^ ^U, 1 . 1 1 , 


c 




\ =(0, 2-2q-r(,, 2) : 


13 


6 




O 

o 






o 


4 




>. 






'w 


2 




c 

o 






c 


n 





2 4 6 8 10 12 
Temperature (K) 

FIG. 3: Temperature dependences of the peak intensities due 
to the Tb ordering at (a) (0,0.424,1) = (0, gTb,l) and (b) 
(0,1.151,2) = (0,2-2gTb,2). 



the peaks at 2qTh remains an open question. 

Next, we measured temperature dependence of a pro- 
file of an ^-type Mn peak to derive the change of the wave 
vector of the Mn spin ordering as a function of tempera- 
ture. The NN-NNN model predicts that the wave number 
shows a step- wise temperature variation as locking at cer- 
tain rational values^ In the previous diffraction measure- 
ments, however, only a smooth temperature dependence 
below Tiock was observedr*^ suggesting that the effect of 
the locking is, if it exists, relatively weak. In the present 
study, we collected the temperature dependence data at 
substantially fine intervals of temperature (AT ~ 0.5 K) 
with the tight coUimation (20'-20'-20'-blank) to detect 
small anomalies in the wave vector of the Mn ordering. 
With this condition, the widths of the Mn peaks arc still 
resolution limited for both the [010] and the [001] direc- 
tions, while the Tb peaks have finite intrinsic widths. 
The correlation lengths of the Tb ordering, which is de- 
rived from the peak width of (0, gxb, 1) peak at 5 K by 
assuming a Gaussian line shape and deconvoluting the 
resolution, is ~140 A both along the [010] and [001] di- 
rections. 

We show in the inset of Fig. ^a) the temperature 
variation of the profile of the Mn peak at (0, gMn,l) ~ 
(0, 0.28, 1). One can see clearly the peak position shifts 
toward higher Q as the temperature increases. By fit- 
ting these data to Gaussians (solid lines in the fig- 
ure) , we extract the temperature dependence of the peak 
height, the peak position (the wave number qun), and 
the peak width (full width at half-maximum, FWHM), 
which are summarized in Fig. ^ With increasing 
temperature, the peak intensity shows an increase at 
T < T^^ ^ 7 K, then decreases [Fig. E^a)]. The 
slope of the temperature profile exhibits an inflection 
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FIG. 4: Temperature dependences of (a) the peak intensity, 
(b) the wave number, and (c) the peak width of the Mn 
(0, (/Mn, 1) peak. Inset of panel (a): Temperature dependence 
of the Mn peak profile at (0, QMn, 1) along the [010] direction. 
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Tiock = 28 K. Finally, the intensity vanishes at 
46 K. For the wave number gMn, it is nearly con- 
stant at 0.276 r.l.u. (reciprocal lattice unit) « 5/18 (or 
11/40) r.l.u. at low temperature with a slow decrease as a 
function of temperature [Fig.^b)]. It shows a steep rise 
above '^34 K, but the upturn of the wave number starts 
at Tiock [inset of Fig. ^b)]. It changes the slope of its 
temperature dependence at 38 K, then it increases mono- 



tonically without any locking behavior. For the peak 
width, it is resolution hmited below T^", and diverges 
above T^" [Fig. EIc)]. We also measured temperature 
dependence of the profile along the [001] direction, but 
did not find any anomalies. 

According to the NN-NNN model, the wave number of 
the Mn ordering gMn decreases in a step- wise manner as 
lowering the temperature, and reaches zero in the ground 
state. However, the step-wise behavior is indistinct and 
qun locks at the finite value below Tiock in the real ma- 
terial, although the observed temperature dependence 
of QMn semiquantitatively agrees with the model. This 
means some higher order effects than the NN and NNN 
spin interactions that stabilize the modulated structure 
must be considered. The coupling between the Mn mo- 
ments and the Tb moments may play a role, but it fails 
to explain the fact that the similar incommensurate spin 
ordering and the locking behavior at the low temperature 
were also observed in orthorhombic YMnOa, where the 
rare earth ion has no magnetic moment >^ Further studies 
are needed to solve this problem, and the study of the 
dynamics will be helpful. 

To sum up the present study, we have investigated the 
ordering of the Mn and Tb moments in a single crystal 
of TbMnOa by neutron diffraction. The ordering pro- 
cess of the magnetic moments are summarized as follows: 
At T^" = 46 K, the Mn moments shows a sinusoidally 
modulated vl-type ordering with a propagation vector 
(0,(7Mn,0), as reported by Quezel et alrL In addition, 
small G, C, and F-type with the same propagation vector 
coexist. gMn is '^0.29 at T^", and smoothly decreases as 
the temperature is lowered showing a weak anomaly at 
38 K. At Tiock = 28 K, gMn locks at 0.276 « 5/18 or 
11/40. The locking of the wave vector amplifies the Mn 
moment, especially enhances the G, C, and F-type com- 
ponents, and squares the spin arrangements. The Tb 
moments order below T'P' = 7 K, which promotes the 
development of the additional Mn components. 
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